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ABSTRACT 

We  consider  the  motion  of  two  rings  of  liquids  with  different  viscosities 
and  densities  lying  between  concentric  cylinders  that  rotate  with  the  same 
angular  velocity  Q.  Gravity  is  neglected  and  interfacial  tension  is  included. 
We  show  that  rigid  motions  are  globally  stable  and  that  the  shape  of  the  In¬ 
terface  which  separates  the  two  fluids  is  determined  by  a  minimizing  problem 
for  a  potential  P  defined  as  the  negative  of  the  sum  of  the  kinetic  energies 
of  two  rigid  motions  plus  the  surface  energy  of  the  interface.  We  shew  that 
the  stable  interface  between  fluids  has  a  constant  radius  when  heavy  fluid  is 
outside  and  J  is  larger  than  one,  where  J  -  -d*[[p]]Q*/T  where  d  is  the  mean 
radius,  [[p]]  <  0  the  density  difference  and  T  the  surface  tension.  When  J  is 
negative  the  heavy  fluid  is  inside  and  the  interface  must  be  corrugated.  The 
potential  of  flows  with  heavy  fluid  outside  is  smaller,  thus  relatively  more 
stable,  than  when  light  fluid  is  outside,  whenever  J  is  large  or  for  any  J 
when  the  volume  ratio  m  of  heavy  to  light  fluid  is  greater  than  one.  These 
results  give  partial  explanation  of  the  stability  and  shape  of  rollers  of  vis¬ 
cous  oils  rotating  in  water  and  the  corrugation  of  the  free  surface  of  films 
coating  rotating  cylinders. 
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SIGNIFICANCE  AND  EXPLANATION 


The  authors  study  the  stability  of  flows  of  two  immiscible  fluids  between 
two  cylinders  rotating  at  the  same  frequency.  They  show  that  this  problem  can 

be  reduced  to  the  minimization  of  a  certain  functional.  Necessary  and  suffi- 

> —  *- 

cient  conditions  for  the  stability  of  "centrifuged"  configurations  with  the 
heavier  fluid  on  the  outside,  are  given.  It  is  shown  that  the  viscosity  plays 
no  role  in  the  stability  problem.  This  is  different  from  shearing  flows 
between  rotating  cylinders,  where  it  is  in  fact  possible  to  have  the  heavier 
fluid  inside. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC  and  not  with  the  authors  of  this  report. 


STABILITY  OF  RIGID  MOTIONS  AND  ROLLERS  IN 
BICOMPONENT  FLOWS  OF  IMMISCIBLE  LIQUIDS 

Daniel  D.  Joseph*'1,  Y.  Renardy2,  M.  Renardy2,  and  K.  Nguyen** 
1.  INTRODUCTION 


We  consider  tne  flow  of  two  immiscible  liquids  with  different  viscosities 
and  densities  lying  between  concentric  cylinders  both  cf  which  rotate  with  the 
same  angular  velocity  ft.  We  neglect  gravity  and  include  interfacial  tension. 
We  study  the  stability  of  steady  rigid-body  rotation  in  which  the  two  fluids 
are  arranged  in  two  rings  with  a  given  volume  ratio.  We  show  that  rigid-body 
rotation  is  globally  stable,  and  the  interface  shape  between  the  two  fluids  is 
determined  by  a  minimizing  problem  for  a  potential  defined  as  the  negative  of 
the  sum  of  the  kinetic  energies  of  two  rigid  motions  plus  the  surface  energy 
of  the  interface.  We  show  that  the  interface  between  the  two  fluids  has  a 
constant  radius  when  the  heavy  fluid  is  outside  and  J  is  larger  than  one, 
where  J  is  defined  in  the  abstract.  This  implies  that  centrifuged  configura¬ 
tions  lie  outside  an  interface  of  constant  radius.  We  note  that  the  stable 
rollers  of  oil  in  water  observed  by  Joseph,  Nguyen  and  Beavers  0  98*0  have 
heavy  fluid  (water)  outside  an  interface  of  essentially  constant  radius.  The 
rollers  are  maintained  in  nearly  rigid  motion  by  the  high  viscosity;  there  is 
no  outer  cylinder,  and  the  motion  of  the  water  is  not  rigid.  We  show  that  the 
Interface  on  rigid  motions  with  heavy  fluid  inside  must  be  corrugated.  Photo¬ 
graphs  of  corrugated  interfaces  of  liquid  films  coating  cylinders  rotating  in 
air  can  be  found  in  the  papers  of  Yih  (I960)  and  Moffatt  (1977).  Yih  (i960) 
gave  a  linear  stability  analysis  for  films  coating  cylinders  rotating  in  air 
and  his  results  are  consistent  with  ours. 

The  potential  of  flows  with  heavy  fluid  outside  is  smaller,  thus  rela¬ 
tively  more  stable,  than  when  light  fluid  is  outside,  when  J  is  large  or  for 
any  J  when  the  volume  ratio  m  of  heavy  to  light  fluid  is  greater  than  one. 
This  is  consistent  with  the  idea  that  configurations  with  heavier  fluid  outside 
should  be  more  stable  because  of  the  centrifugal  force  and  that  if  the  inner 
fluid  is  heavier,  the  rigid  motion  should  be  less  stable  no  matter  what  the 

*  Department  of  Aerospace  Engineering  and  Mechanics,  University  of  Minnesota, 
Minneapolis,  MN  55455. 

**  Firestone  Tire  and  Rubber  Co.,  Akron,  OH  44317. 


^ponsored  by:  ’  ~~ 

The  Fluid  Mechanics  Branch  of  the  National  Science  Foundation  and  the 


United  States  Army  under  Contract  No.  DAAG-29-82-K0051 . 

The  United  States  Army  under  Contract  No.  DAAG29-80-C-004 1 .  This  material 
is  based  upon  work  supported  by  the  National  Science  Foundation  under  Grant 
No.  MCS-82 10950  and  8215064. 


b  q? ^ wwwwwwwy'P ]■. -,  •)  WMIHWH  ,*  sw  ■■  i.y. rr-r  r-rr 


8 


-2- 

viscosities.  The  stability  of  flows  in  which  viscosity  differences  are  impor¬ 
tant  depend  strongly  on  the  viscosity  ratio.  Rigid-body  rotation  involves  no 
shear.  As  a  result,  the  mechanism  which  we  call  lubrication  stabilization 
(Renardy  &  Joseph,  198M),  in  which  thin  layers  of  the  less  viscous  fluid  occupy 
regions  of  high  shear,  is  absent. 

2.  EQUATIONS  OF  MOTION  AND  INTERFACE  CONDITIONS 

Consider  the  flow  of  two  immiscible  liquids  contained  between  two  infin¬ 
ite  concentric  cylinders.  The  perturbed  regions  occupied  by  liquid  1  and  2  are 
denoted  respectively  as 

V,(t)  -  (r.e.xja  i  r  i  R(x,8,t),  -•  <  x  <  •,  OSes  2*), 

V,(t)  -  (r,0,x|R(x,8,t)  SrSb,  -•  <  x  <  •,  OSes  2»}. 

The  stress  is  given  by 


T  -  -pi  ♦  S,  S  -  2U  D  [u ]. 


The  equations  of  motion  hold  in  each  region 


(2.1) 


div  U£  -  0,  l  -  1,2, 


(2.2) 


Pi  dt  Hi  "  "V*  4  div  £l» 


(2.3) 


where  u£  -  ej.u£  ♦  e0V£  ♦  exW£,  ♦  -  p  ♦  pgr  sine,  D[u3  -  ^r(?u  ♦  VuT),  p,  and  pa 
are  the  densities  and  ui  and  ya  are  the  viscosities.  In  all  that  follows  g  -  0. 
The  cylinders  at  r  -  a  and  r  -  b  rotate  with  some  constant  angular  frequency 
G.  At  the  interface  Z 


F(x,r,e,t)  -  r  -  R(x,6,t)  -  0, 


(2.*0 


j% 

f! 

1 
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we  have 


dR 


v  3R 


3R 
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(2.5) 
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We  also  require  that  the  Jumps  across  £ 


[M3  -  (•),  -  (•), 

In  the  velocity,  the  shear  stress  and  the  difference  between  the  jump  in  the 
normal  stress  and  the  surface  tension  force  all  vanish. 

We  are  going  to  study  spatially  periodic  solutions  which  are  2*/a  periodic 
in  x  and,  of  course,  2n  periodic  in  e.  The  volume  of  each  component  fluid  is 
prescribed  by  specifying  a  mean  radius 

d*  -  R*  (2.6) 


Our  convention  is  the  fluid  with  subscript  1  is  on  the  inside.  The  fluid  on 
the  Inside  can  be  heavy  or  light.  We  are  Interested  in  two  cases 

(A)  The  heavy  fluid  inside,  a  i  r  S  R*(8,x,t) 

CCp33  >  0,  d  *  -  R*  (2.8) 

(B)  The  heavy  fluid  is  outside,  Rg(8,x,t)  i  r  $  b 

CCp33  <  0,  d*  -  R*.  (2.9) 


The  volume  of  light  fluid  and  the  volume  of  heavy  fluid  is  fixed,  indepenaent 
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of  whether  it  is  inside  or  outside  when  mA  -  mg.  Then 


2  2 

a  m.+a*  .  b*  ♦  in 

M  _ b _  M  _ a 

dA  1  ♦  m  ’  dB  1  ♦  m 


(2.11) 


The  Jump  equations  at  r  •  R  are 


CCu]]  -  0, 


“(CCpD]  ♦  2HT)n  ♦  [[§]]  •  n  -  0, 


(2.12) 


where 


RRee(TR*)  ♦  RRXX(R**R*)  -  R*(1*R*)  -  2R*  -  2RR0RXRX0 
(R*  ♦  R  ♦  R2R  )%/a 

O  X 


(2.13) 


2H  is  the  sum  of  the  principal  curvatures,  T  is  the  surface  tension,  and  n 
VF/|VF|  and  VF  -  £p  -  e0R0  -  2x^xm 

3.  RIGID  ROTATION  OF  TWO  FLUIDS 


The  velocity 


u,  -  e0Qr, 


(3.D 


and  the  pressure 


(3.2) 


where  the  constants  (p,D)  are  (p|PD,)  in  V,  and  (pa,Da)  in  Va,  is  a  solution  of 
(2.2),  (2.3)  and  (2.4)  with  S  identically  zero.  We  suppose  that  R(8,x,t),  peri¬ 
odic  in  x  and  8,  is  prescribed  and  arbitrary.  At  r-R  we  have 


v  v  V  .'  v  .■ 


“/  V  V  V  r .  r ,  ,  .  .  .  .  . 


*,  s 


Ctue]]  -  0, 


Up «]j  •  Up]]  f  R*  ♦  [[D)J. 


(3.3) 


We  cannot  satisfy  the  differential  equation 


CCp#33  ♦  2HT  -  0, 


(3.*0 


expressing  the  Jump  condition  for  the  normal  component  of  the  stress,  for  an 
arbitrary  given  surface  R(6,x,t).  We  call  (3.1)  ”  (3*3)  an  "extended"  rigid 
motion  and  we  prove  that  these  motions  are  globally  stable  with  the  shape  of 
the  interface,  and  possibly  the  placement  of  heavy  and  light  liquid  determined 
by  a  new  minimum  principle  of  classical  type. 

4.  PERTURBATION  EQUATIONS 

Let  u,p  be  the  velocity  and  pressure  in  the  deformed  domain  and  let 

u  -  u,  ♦  u, 

P  •  Po  ♦  P 

M  m, 

where  u0  -  e6V(r),  and  pc(r)  is  the  associated  pressure.  The  function  u,  p  and 

A 

S  are  perturbations.  In  designating  components 

u  »  (u,v,w) 

we  suppress  the  caret  overbar.  All  these  quantities  are  defined  in  V,(t)  and 
V,(t).  For  the  moment  we  leave  open  the  possibility  that  [[V]]  i  0,  (for  rigid 
motion  [[V]3  -  0).  The  equations  governing  u  and  p  are 


3u  a  .  .  .  .  . 

p[(^  *  -  *  ^-°  +  “  ~^P  ♦  V  •  S, 


(«U) 


V  •  u  -  0 


(^.2) 


J 


The  boundary  conditions  are 


u,(r-a)  -  u,(r-b)  -  0. 


(U.3a,t>) 


The  interface  conditions  on  r  •  R(x,6,t)  are 


CCu]]  -  [[«]]  -  CCvV]]  -  0 


(!».»») 


9R  v+V  9R 

U  *  9t  R  96 


(*».5) 


-CCp]]n  ♦  CCS]]  •  n  -  CCp.]]n  ♦  2HTn 


<*.6> 


5.  ENERGY  EQUATIONS  FOR  NONLINEAR  DISTURBANCES 
Introduce  the  following  notations: 


<•>  •  J  (0  dV, 


-  |  (*)d  l. 


where 


dt  -  RdedxUUR  */R*+R* 


Rd6dx| VF| 


and  we  are  assuming  that  the  disturbance  flow  is  2*/a  periodic  in  x. 

To  form  the  energy  equation  we  multiply  (*».  1 )  by  u,  integrate  over  V,  and 
V„  add  and  use  Reynolds  transport  theorem  to  show  that 


I 


*•.  *\ 


V*, 


and 


♦  ^pu  •  KG*]  •  u)  ♦  ^2yD[u]:D[u]  ^ 


^C[-pH‘  •  n  ♦  u  •  S  •  n]]^£ 


(5.1) 


For  flow  in  circles  with  u0  -  e^rO  we  have 


^pu  •  DCu0]  •  u^-  ^puvr 


which  vanishes  on  rigid  motions.  Moreover,  since  u  •  u  -  u0  is  continuous 
across  I,  we  have 


^[(u  •  (-pn  ♦  S  •  n)]]^E  -  ^u-u0)  •  [[-pn  ♦  S  •  n]]  ^ 
-  ^(u-Uo)  •  n([[p,]]  ♦  2HTj^£ 


(5.2) 


where  the  last  equality  follows  from  (4.6).  We  next  observe  (see  Joseph,  II, 
Eq.  (96.11))  that 


(*!  •  i)l  ■  '5 


TdE. 


(5.3) 


Since  T  is  constant,  we  have 


\  * 

Li 

& 

•  . 

<m  ' 


■c** 

k: 

IC/' 

tv 

Kv 


1/  *  ■ 


■s 


r* 

>.* 

V  ■- 


>\ 

k.% 

*„\ 


&7- 


£ 

v* 

V’ 

Cv 

<V 

»V 

& 

P\  . 

» v%.- v  . 
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^2HTu  •  n  ^  -  -T  ^((^Ka«»K  *  ♦  R**  *  )) 


where  ((•))  is  defined  in  (2.7). 


Moreover,  using  (2.5)*  we  find  that 


VF  1  3R 

u  •  n  ■  u  •  - — -  -  : — 7  sr. 
“  ~  _  VF  VF  3t 


Hence 


^CCp,]]u  •  -  «[[p0]]R  |5)). 


Since  CCp#]]  -  CCp]]QaRa/2  ♦  CCD]]  is  a  function  of  R  alone,  there  is  a  scalar 
function  ♦(R), 


♦(R)  -  CCp]]Oa  --Y— » 


(5.H) 


such  that 


((CCpo]]R  §|))  -  ~  ((♦(R))); 


in  deriving  (5.*0  -  (5.5),  we  used  (2.6)  to  set 


at  ((fl'»  '  °- 


(5.5) 


Finally,  we  show  that,  for  rigid  motions  where  ue  •  n  -  &Rg/|VF|,  we  have 

*  £(C[po33  ♦  2HT)^j. 


-  n((RRe(CCp03]  ♦  2HT)))  -  0 


(5.6) 
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Ef 


In  deriving  the  first  equality,  we  use  n  -  7F/| VF | ;  then  we  note  that  since 
[Lp0]]  is  a  function  of  R  alone,  there  is  a  2*  periodic  m(R)  such  that 

MCCp#})  -  3m(R)/ae 


which  vanishes  on  integration.  We  next  use  the  expression  (2.13)  to  write 


£ 


k.v 


Since  R(e,x,t)  is  periodic  in  x  and  0,  the  last  integral  in  (5.6)  vanishes. 
Collecting  all  these  results,  we  find  that 


(5.7) 


where 


E  - 


H 


D  -  ^uDCu]:D[u]), 

T  *  *1  ,  [[p]]0* 

T((|r**R  g  *  R*R  XJ1^*)) - s —  ((CR#"d*]*)). 


(5.6) 


The  function  H[R]  is  the  variable  part  of  a  "potential  energy"  P[Rj  for  rigid 
motions.  It  is  easily  verified  that 


1  .  *  ,  ■  .  •  .*•*.*/  /  «•  *  ’  «  •  '  •  •  V  ’  k  •»  *  w  "  ■  *  _  ■  _*  %  .  %  .  V  _%  .  V  %  *,  *,  O  O  W  •  „ 


'  k 


V"  ^  v.  V"  -.  V^v  'S.  V*  % 


jMJ 


K»  -  [-((a*))  ♦  ((R^Jjp.fiVS, 

Ka  -  [((b*))  -  ((R^np.QVs 

is  the  kinetic  energy  of  a  rigid  notion  in  region  V,  and  Va,  respectively.  The 
associated  potential  of  this  rigid  notion  is  the  negative  of  the  kinetic  ener¬ 
gies  plus  the  surface  energy 


P  d®f  -(Kj*Ka)  |  dE 
c»  ♦  T(([Rao*R*)*R*],/a»  -  tcpttdR'W/e 


C*  ♦  M 


(5.9) 


where  Ct  and  Ca  are  different  constants.  We  nay  write  (5.7)  as 


d(E*P) 

dt 


(5.10) 


Equations  (5.7)  and  (5.10)  were  derived  by  D.  D.  Joseph. 

In  the  next  section  we  shall  show  that  rigid  notions  are  globally  stable, 
as  is  the  case  with  one  fluid,  but  that  the  stable  configurations  of  the  rigid 
notions  nininize  P  subject  to  the  volune  constraint  (2.6). 

It  is  useful  to  write  the  potential  in  a  dimensionless  forn  in  which  R  ■>  d 
♦  6  and  (6,x,1/a)  are  made  dimensionless  with  d  and  A  -  6/d.  Then,  to  within 
constants,  we  may  define  a  dimensionless  potential 

P  -  -((D)  ♦  (([(1*A)a(1*A*)  ♦  A*],/a))  ♦  jj(([2A* A*]a)),  (5.11) 


where 


.YJ  ‘  iOvIy  1  v 
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,  d*[[p33o* 

j  - — - 

has  the  sign  of  CC p 3D.  The  constraint  (2.6)  implies  that 

<(2A  ♦  A1))  -  0.  (5.12) 

This  shows  that  the  average  deviation  A  -  -1  ♦  R/d  from  zero  must  be  negative 
if  the  volume  of  the  two  fluids  is  preserved. 

It  Is  necessary  to  remark  that  the  representation  r  -  R(e,x,t)  of  the  free 
surface  is  not  completely  general  and  it  loses  Its  utility  when  the  magni¬ 
tude  |  &  |  of  the  deviation  from  the  cylinder  is  equal  to 

A  -  minjl  -  |,  |  -  l].  (5.13) 

If  )a|  2  A,  then  the  interface  will  touch  one  or  the  other  of  the  cylinders, 
and  a  smooth  free  surface  will  not  be  possible. 

The  linearized  form  of  (5.11)  for  A  near  to  zero  is 


P  -  |  ((CJ-13A*  ♦  A*  ♦  A*)). 


(5.1^) 


6.  STABILITY  OF  RIGID  MOTION 

The  following  results  concern  solutions  of  the  equations  which  are  smooth 
for  all  time. 

Theorem  1 .  Rigid  motions  are  stable  in  the  sense  that  periodic  distur¬ 
bance  of  rigid  motion  must  decay  in  the  mean.  We  first  note  that  D  >  0  on  all 
non  zero  disturbances  of  rigid  motion.  We  are  considering  the  stability  of 
flows  with  heavy  fluid  outside,  -[[p]]  -  pa  “  Pi  >  0,  or  inside,  [Cp33  >  0. 
These  two  situations,  called  (A)  and  (B)  in  Section  2,  are  distinct  in  that  they 
cannot  be  connected  by  time  dependent  motions  with  smooth  interfaces.  We 
choose  u  to  be  a  disturbance  of  one  or  the  other  of  these  two  possibilities. 
Then  the  positivity  of  D  Implies  the  decay  of  E  ♦  P.  In  fact,  we  could  show 
that  there  is  X  depending  on  the  viscosities  and  the  densities,  such  that  D  i 
IE  where  X  depends  on  v„  v„  p„  p„  y,,  y,.  Integrating  (5.10)  from  t  •  0  to  t, 
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we  find  that 


E(t)  ♦  P(t)  -  E(0)  ♦  P(0)  -  D( t)  dt 

Jo 


r 

i  E0  ♦  P,  -  Jtj  E(t)  dt;  (6.1) 


it  follows  that 


f 

1  E( t)  d-t  S  ECO)  ♦  P(0)  -  E(t)  -  P(t).  (6.2) 

Jo 


Since  P  is  bounded  rrom  below  we  conclude  that 


E(t)  and  D(t)  are  integrable.  (6.3) 

Moreover,  assuming  that  E  goes  to  0  as  t  ■*  ♦  •,  P  admits  a  finite  limit  as  t 
goes  to 

Let  us  consider  the  limit  configuration  (£(•), P(«))  since  E(»)  -  0,  this  is 
a  rigid  motion.  To  show  that  P(«)  is  a  minimum  of  the  functional  P  as  R 
varies  we  consider  any  rigid  motion  (E(0)  -  0,  P(0) )  and  assume  that  this  rigid 
motion  goes  to  the  rigid  motion  (E(«)  -  0,P(»))  as  t  goes  to  ♦«*.  If  P(0)  i 
P(«),  then  the  interface  between  the  two  liquids  must  have  moved  from  the  con¬ 
figuration  at  t  -  0  to  the  one  at  t  ■  •,  and  in  this  motion  D(t)  >  0  on  an  in¬ 
terval  with  non  zero  measure.  Then,  from  (5.10) 


P(»)  -  P(0)  -  -  D(t)  dt  <  0,  (6.H) 

Jo 


so  that 


P(-)  <  P(0) 


(6.5) 


Thus  P  decreases  ir  every  change  of  configuration  between  rigid  motions.  Since 
P  is  a  bounded  from  below  functional  of  R,  P  must  decrease  to 


p<-)  •  Jiipu>  -  BRinp<*>. 


(6.6) 


where  R*  •  d*  and  R  is  periodic  and  continuously  differentiable  in  x  and  6. 

We  may  describe  the  set  of  such  periodic  functions  R  with  E*  •  d*  as  a  one 
parameter  family  R(c).  We  assume  that  R(0)  -  R  where 


“Rn  P(R)  -  "cln  P(R(c)). 


(6.7) 


We  must  have: 


dl  CR(c)]  c-o  *  "((R  3I(0)(2RT  *  -  0. 


(6.8) 


where  2H  denotes  the  curvature  (2.13)  evaluated  on  the  minimizer  R  -  R.  Since 
l*(c)  -  da,  we  have 


((R  57(0)))  -  0, 


(6.9) 


»  dR 


so  that  R  37(0)  Is  orthogonal  to  constants.  It  follows  now  from  (6.8)  and 
(6.9)  that 


2HT  ♦  CCpJ]  -  C  -  2HT  ♦  [[p]]  (6.10) 

where  H  is  the  mean  value  of  H. 

Equation  (6.10)  may  be  recognized  as  the  differential  equation  arising 
from  the  normal  stress  condition  (2.12)  on  rigid  motions. 

It  now  follows  that  extended  rigid  motions,  which  are  globally  stable, 
are  actually  hydrodynamically  admissible,  with  a  balanced  normal  stress  equa- 
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tion,  when  the  Interface  R  -  R  is  a  minimizer  of  (6.6).  The  stable  rigid 
motions  satisfy  (2.5)  in  the  form 


0 


3R 

36* 


where  R  is  time  independent  in  a  rotating  coordinate  system. 

Theorem  2.  The  stable  configurations  are  those  rigid  motions  which  min¬ 
imize  P  among  C'(x,8)  functions  R(x,6)  satisfying  the  volume  constraint  (2.6). 

Consider  first  case  (A)  in  which  the  heavy  fluid  is  outside,  J  >  0.  If  J  > 
1,  then  P  given  by  (5.14)  is  a  minimum  when  A  -  0. 

Theorem  3.  The  cylindrical  Interface  with  constant  radius  R  -  d,  is 
stable  against  small  disturbances  if  and  only  if  J  2  1. 

It  is  of  interest  to  ask  when  R  -  d  is  a  global  minimum  of  P  among  all 
interfaces  of  the  form  r  •  R(e,z)  compatible  with  the  volume  constraint.  This 
question  is  answered  by  the  following  theorem  of  M.  Renardy. 

Theorem  4.  The  concentric  Interface  R  -  d  is  a  global  minimum  of  P  among 


all  interfaces  r 


R(e,z)  satisfying  (3.4)  if  and  only  if  J  2 


First  we  show  that  the  criterion  of  theorem  4  is  sufficient  for  stability. 


Certainly  we  have  P2P-  TR  ♦  3  O^Pa-p^R*  ,  and  it  is  sufficient  to  show 
that  R  -  d  minimizes  P.  We  set  R*  -  d*(1  +  Y),  hence  Y  is  subject  to  the  con- 
straints  *1  ♦  jj  S  Y  S  and  ((Y))  •  0.  We  then  have 


y 


We  define 


P  -  Td  ((/ 1  ♦Y  ♦  j  J  (1*Y)*)) 

•  Td  ((/UY  ♦  |  J  (1*Y)*  -  y[1  ♦  J  j])).  (6.11) 

•  i J  «'•«*  -  t[j  *  t  ■>]• 


f ( Y)  •  /I  ♦  Y 


(6.12) 


The  constant  multiplying  Y  has  been  chosen  such  that  f'(0)  -  0.  If  J  2 

<4 

— - rj,  then  f  has  its  minimum  at  T  •  0,  for  Y  in  the  range 


[<  *  i] 

[-1  ♦  |5.  *1  ♦  §»]- 


The  criterion  is  also  necessary  for  a  global  minimum.  He  can  choose  R 
such  that  Rg  -  0,  and,  by  choosing  a  small  (long  waves)  we  can  make  Rz  as 
small  as  we  like.  It  follows  that  it  is  also  necessary  that  R  -  d  minimizes 


P.  If  J  < 


[■  •  a 


r  i  then  T 


[-  •  ad 


<  f (0),  and  the  graph  of  f  is  sketched  as 


a*  a* 

Figure  t.  If  we  draw  the  tangent  from  the  point  (-1  ♦  f(-l  ♦  ^j))  as  indi¬ 
cated,  it  will  touch  the  graph  at  a  point  a,f(a))  to  the  right  of  Y  -  0.  and  it 
will  intersect  the  line  Y  -  0  at  a  value  below  f(0).  This  means  that  we  can 

a* 

take  a  perturbation  Y,  alternating  between  the  values  -1  *  and  a,  such  that 
(( Y))  -  0.  but  ((f(Y)))  <  f(0). 

Turning  next  to  case  B  in  which  the  heavy  fluid  is  inside  J<0  we  find  that 

•  a  2 

f ( Y)  is  concave.  Hence  P  will  have  minimizers  only  at  boundary  values  -1  ♦  ^ 

and  -1  ♦  In  an  infinite  cylinder,  there  will  therefore  be  no  minimizers  of 

P  of  the  form  r  -  R(8,z).  In  a  finite  cylinder,  we  cannot  make  Rz  arbitrarily 
small  without  also  making  Y  small,  and  there  can  be  stable  motions  with  heavy 
fluid  inside,  which  have  a  corrugated  free  surface,  as  in  the  experiments  of 
Ylh  (I960)  and  Moffatt  0977).  It  would  be  of  interest  to  determine  these  cor¬ 
rugated  shapes  as  a  solution  for  the  minimum  problem  of  P.  If  J  is  large,  we 
expect  the  amplitude  of  the  corrugated  surface  to  be  also  large  and  eventually 

a*  ba 

violate  the  constraints  -1  +  ^s<Y  <  -1  ♦ 

The  results  of  this  paper  have  some  relevance  for  the  problem  of  centri¬ 
fuging.  Intuitively  one  expects  that  the  heavy  fluid  will  be  outside  if  the 
rate  Q  of  rotation  is  large,  even  if  the  heavy  fluid  were  initially  on  the 
inner  cylinder.  The  transport  of  fluid  from  the  inner  to  the  outer  cylinder  is 
a  topologically  complex  process  which  cannot  be  handled  in  the  frame  of  our 
smooth  parameterization  of  the  interface.  The  transport  of  fluid  from  the 
inner  to  the  outer  surface  is  also  a  physically  complex  process  involving  the 


\V\- 


.  ••  •• 
.'A-*.--/.*.'- 
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rupture  of  adhesion  at  both  walls  and  possibly  internal  fracturing  and  healing 
of  the  liquids  themselves.  These  physical  processes  are  not  well  understood, 
and  they  do  not  appear  in  our  equations.  Nevertheless,  it  is  not  unreasonable 
to  seek  stable  configurations  among  those  which  minimize  P  with  respect  also 
to  the  position  of  heavy  and  light  fluid. 

He  can  compare  the  potentials  P A  and  Pg  for  flow  with  heavy  fluid  inside 
t[p]]  >  0  and  heavy  fluid  outside  [Cp33  <  0  under  the  condition  that  the 
volume  m,  defined  in  Section  2,  is  fixed.  The  total  kinetic  energy  KA  wherv  the 
heavy  fluid  with  density  pt  is  inside  is  given  by 

p.Q*  %  p,Q*  * 

KA  *  Kj  ♦  Kj  -  -  g—  ((R^-a*))  *  — g —  (Ob'-RR^)) 

a  a 

where  ((RA))  •  (<d A )).  When  the  heavy  fluid  with  density  pt  is  outside 

p.Q*  *  p.Q*  * 

Kg  ■  ♦  Ka  ■  — g—  ((b*~Rg))  ♦  —g—  ((Rg“a,‘)) 


where  (CR B ) )  •  ((dg)).  Using  (2.10),  with  mA  •  mg  ■  m  we  find  that 


The  potentials  are 


where 


PA  *  “kA  *t  |£a|» 
PB  *  -Kg  *T  | Eg | 


■I.: 


. 

%  «• « a  xoT- ^4  A/,  ■>,  f.  «»„  n'.v,,  •_  •*,  •*.  a  -•  **  »*  .»  .■  ,»  v  .* 

V  v  V  v  V  r  j*  *>  v  '  »  V  v  *-•  V  *.•  *«*%.*\k  *  -  *  *  ^  .**  .*•  «,*»  .’** 


vv 


Is  the  area  of  IA.  The  difference  in  the  potentials  is 

Pa  -  Pb  -  -ka  ♦  kb  ♦  t(|ea|-|ib|) 

fCflTJQ*  *  * 

-  <(b**a*-RA-RB)>  ♦  T(|rA|-|EB|) 

Here  [[p]]  -  p,  -  pa  >  0.  If  PA  >  PB  then  (B)  with  heavy  fluid  outside  is  more 
stable.  Since  ((R2))  -  ((d2)),  we  have 

(([R*-da]a))  -  <(R*))  -  ((d*)). 

The  potential  difference  may  be  written  as 

'k-  rB-  ((u-'i'-aj-aj))  -  *  [v^b}”]  *  ^M'M* 

*  "" 


(6.1) 


Now,  using  (2.11),  we  eliminate  dA  and  d0 


Pa  -  Pb 


-  (<[Ri-s]’  •  [«;-<]■»]  *  tie*i  -  tm 


(6.2) 


Consider  now  the  case  of  uncorrugated  interfaces  with  RA  -  dA,  RB  -  dB,  EA 
2 idA,  Eb  -  2vdg.  We  find  that 


Pa  -  pb 


CCp]]Q*  (([b^a2]2)) 


2wT((dA-dB)) 


(6.3) 


where  dA  and  dB  are  given  in  terms  of  m  by  (2.11).  If  the  volume  ratio  m  of 
heavy  ( p ,)  to  light  (pa)  fluid  is  greater  than  one,  then  dA  >  dB.  Hence  PA  >  PB 
when  m  >  t.  If  there  is  only  a  small  amount  of  heavy  fluid,  m  <  1  and  dB  <  dA 
then  PA  >  PB  if  CCp]]Q*/T  is  large.  In  all  these  cases  the  configuration  with 
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where  dA  and  aB  are  given  in  terms  of  m  by  (2.11).  If  the  volume  ratio  m  of 

heavy  (p,)  to  light  (p2)  fluid  is  greater  than  one,  then  dA  >  dg.  Hence  PA  >  Pg 

when  m  >  1.  If  there  is  only  a  small  amount  of  heavy  fluid,  m  <  1  and  dg  <  dA 

then  PA  >  Pg  if  C  [  p 3 3Qa/T  is  large.  In  all  these  cases  the  configuration  with 

heavy  fluid  outside  is  more  stable.  If  [[p]]Q*/T  is  small  enough,  then  PA  < 
Pg,  and  the  configuration  with  heavy  fluid  inside  is  more  stable.  However  our 
earlier  analysis  showed  that  when  the  heavy  fluid  is  inside,  the  Rg(e,x)  which 
minimizes  P  is  not  everywhere  equal  to  dg. 

7.  STABILITY  OF  ROLLERS 

Rollers  are  viscous  fluid  bodies  which  rotate  as  rigid  wheels  in  fluids  of 
smaller  viscosity.  These  rollers  have  been  observed  (Joseph,  Nguyen  and 
Beavers,  1983)  in  bicomponent  flows  of  immiscible  liquids  in  several  different 
flow  configurations:  on  a  cylinder  rotating  in  a  box,  between  the  four  cyl<* 
inders  of  Taylors  mill  for  studying  straining  flows,  and  separating  dynamically 
driven  Taylor  vortices  between  rotating  cylinders. 

The  most  interesting  feature  of  the  dynamics  leading  to  the  formation  of 
rollers  is  the  fracturing  of  the  viscous  liquid  at  some  critical  level  of  the 
stress.  In  this  process  the  roller  breaks  away  form  the  side  wall  and 
relieves  the  high  stress  associated  with  no  slip  at  the  side  wall.  So  in  the 
final,  stable  dynamics,  rollers  are  lubricated  by  water  and  air  on  all  sides. 
The  rollers  rotate  nearly  as  rigid  bodies  because  they  are  so  viscous.  The 
stability  of  rollers,  as  our  analysis  suggests,  depends  on  the  fact  that  the 
density  stratification  is  such  as  to  prevent  the  centrifuging  of  the  roller. 
The  viscosity  ratio  is  probably  not  an  important  factor  in  the  dynamics  of 
stable  rollers. 

The  water  which  surrounds  the  rollers  in  experiments  is  at  rest  near  the 
tank  wall  and  cannot  rotate  rigidly.  Therefore  rollers  are  not  a  special  case 
of  rigid  motions  studied  in  this  paper.  However  the  density  stratification, 
with  water  outside,  does  contribute  to  the  stability  of  rollers,  with  a  stabile 
izing  term  C C pV*3 36*  at  the  interface,  where  V  is  the  common  velocity  of  fluid 
particles  on  either  side  of  the  interface  and  6  is  the  surface  deflection;  here 
assumed  small. 

Rollers  are  unstable  to  non-ax isymmetric  disturbances  when  the  angular 
velocity  is  high  enough.  This  instability  is  associated  with  viscous  shearing, 


which  becomes  important  at  higher  speeds  and  with  a  possible  unstable  distri¬ 
bution  of  angular  momentum. 

The  low  speed  rollers  are  robustly  stable.  In  our  analysis  we  did  not 
consider  gravity,  but  gravity  does  enter  into  the  dynamics  of  the  stable  rollr 
ers  reported  in  the  paper  of  Joseph,  Nguyen  and  Beavers  ( 1 9814)  and  here.  In 
experiments  in  which  the  top  of  the  roller  rotates  in  air  the  roller  would 
centrifuge  out  into  the  air  were  it  not  for  gravity  wMch  on  the  small  top 
portion  of  the  roller  exposed  to  air,  is  nearly  radial.  A  similar,  but  smaller 
effect  due  to  gravity  occurs  at  the  bottom  of  the  roller  which  is  pushed  up  by 
gravity  because  the  lighter  oil  is  buoyant  in  water.  Gravity  tends  to  flatten 
rollers  into  right  circular  cylinders.  To  a  degree  the  diameter  of  stable 
rollers  can  be  controlled  by  gravity,  with  a  tendency  fo  the  roller  to  poke 
its  head  into  the  air.  We  are  able  to  change  the  diameter  of  the  rollers  by 
changing  the  water  level  in  the  box.  This  effect  of  gravity  is  exhibited  in 
plates  7.1  a«c.  Sketches  of  the  side  view  of  these  plates  is  shown  in  Figure 

7.1  a»-c. 

The  principal  effect  of  gravity  may  be  eliminated  by  submerging  the 
roller  entirely  in  water,  as  in  Plates  7.2(a)  and  (b).  When  the  flattening 
effects  of  gravity  are  absent  the  shape  of  the  interface  on  stable  rollers  is 
strongly  influenced  by  interfacial  tension,  with  bounding  surfaces  in  nearly 
circular  arcs,  as  in  Figure  7.2(b).  The  pressure  distribution  in  the  water  is 
not  a  strong  barrier  to  centrifuging  and  the  dynamics  of  the  roller  in  Plate 

7.2  are  closer  to  case  B  of  this  paper  in  which  the  heavy  fluid  is  on  the 
inside  cylinder  with  a  corrugated  interface  separating  the  two  liquids. 
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Plate 

7.1(a) 


Plate  7.1:  Holler  of  silicone  oil  (p  -  0.95  gm/cc,  p  -  95,000  cp)  in  water  at 
different  water  levels.  The  rod  is  made  of  plexiglass,  2  inches  in  diameter, 
rotates  at  10  RPM. 


(a)  The  roller  is  very  nearly  in  a  solid  body  rotation  with  small  shearing  by 
water  at  the  roller  rim.  Part  of  the  roller  is  in  water  and  the  other  in  air. 
The  roller  is  very  stable,  held  together  by  hydrostatic  pressure  in  water  and 
gravity  in  air. 
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Plate 


(b)  Water  is  added  to  the  box.  The  roller  becomes  larger  by  flattening  out 
but  remains  round  and  stable. 
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Plate 

7.1(c) 


(c)  More  water  is  added.  The  roller  becomes  even  larger.  The  roller  is  now 
completely  submerged  in  water  and  is  slightly  out  of  round  due  to  buoyancy. 
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angular  velocity  Q.  Gravity  is  neglected  and  interfacial  tension  is  included. 
We  show  that  rigid  motions  are  globally  stable  and  that  the  shape  of  the  in- 
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20.  ABSTRACT  (cont.) 


terface  which  separates  the  two  fluids  is  determined  ly  a  minimizing  problem 
for  a  potential  P  defined  as  the  negative  of  the  sum  of  the  kinetic  energies 
of  two  rigid  motions  plus  the  surface  energy  of  the  interface.  We  shew  that 
the  stable  interface  between  fluids  has  a  constant  radius  when  heavy  fluid  is 
outside  and  J  is  larger  than  one,  where  J  -  **d*[[p]]Ql/T  where  d  is  the  mean 
radius,  [[p]]  <  0  the  density  difference  and  T  the  surface  tension.  When  J  is 
negative  the  heavy  fluid  is  inside  and  the  interface  must  be  corrugated.  The 
potential  of  flows  with  heavy  fluid  outside  is  smaller,  thus  relatively  more 
stable,  than  when  light,  fluid  is  outside,  whenever  J  is  large  or  for  any  J 
when  the  volume  ratio  m  of  heavy  to  light  fluid  is  greater  than  one.  These 
results  give  partial  explanation  of  the  stability  and  shape  of  rollers  of  vis¬ 
cous  oils  rotating  in  water  and  the  corrugation  of  the  free  surface  of  films 
coating  rotating  cylinders. 


